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ABSTRACT 


The broad aim of this thesis is to inquire into the 
nature of the Fourier series representation for shape. 
Specially, the emphasis is upon arriving at a suitable basis 
for classifying shapes which are closely related; the ease 
with which we can accomplish this task is an indication of the 
ability of the representation to capture fine details of the 
shape and to dissociate them from the less varying or gross 
features of the shape. We used two methods to classify 
shapes. In the first, the concept of bending energy stored in 
an elastic rod which is deformed to the shape in question is 
used. The expression for bending energy corresponds to the 
integral for the curvature over the entire boundary and is 
therefore intuitively satisfying. The second method treats 
the descriptors as parameters in a filter and shapes are 
classified by matching the parameters of the target with the 
known parameters as is done for matched filters in signal 
processing. 



CHAPTER 1 


INTRODUCTION 

One of the main areas of interest in artificial 
intelligence is the use of representations for problem 
solving. A representation for an object is a symbolic 
abstraction of the object in the real world. In applied 
mathematics, for example, numbers and operations on numbers 
are used to represent the physical world and the actions in 
it, The main reason why one chooses to deal with represen- 
tations rather than with the actual objects these repre- 
sentations stand for is that they bring important information 
to the foreground, and hence make it easy for us to solve 
the problem in the representation world. The nature of the 
problem dictates the choice of a particular representation; 
in general, the same object can have different representations 
emphasising different aspects of the object in the real 
world. For example, the so-called intrinsic representation 
(curvature vs, arc-length) for the boundary of an object 
brings out the curvedness or wiggliness of the boundary. The 
medial axis transform representation, on the other hand, 
brings out the correspondence of the opposite sides of a 
boundary, thus emphasizing width perception relating boundary 
arcs all along parallel sides. 
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The existence of many representations for an object 
immediately poses many important questions such as the 
extent to which a representation corresponds with the 
real v/orld, of how powerful the representation is and what 
class of objects are representable* 

It was this desire to probe into the nature of the 
representation that led to the present work. To be more 
specific, a pattern recognition system with Fourier descri- 
ptors being used as a representation for a 2D-’Shape had 
been built by Karnick [4 ] , The present thesis was conceived 
of as being complementary to the above work in as much as 
we decided to investigate the adequacy of the Fourier series 
representation for describing shapes-specif ically those of 
hand written characters - in more detail, 

1.1 STATElvEiTT OF THE PROBLEM: 

Conventional pattern recognition is concerned more 
with classification among objects which are generally quite 
different from each other than with the particulars of 
differentiation between objects whose general shape is 
similar. In the recognition of hand written characters, 
for example, the problem of reliably distinguishing between 
two characters, say, an A and B has, understandably perhaps. 
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received much more attention than that of being able to 
distinguish between the natural variants of a character. The 
emphasis has generally been on interclass classification 
rathex than intraclass separation. 

The main concern of this thesis will be with the 
problem of distinguishing between 2D-shapes (we will clarify 
soon what is meant by shape) which, though they resemble 
each other in a global (overall) sense are neverthless, 
recognized by us as being different from each other because 
there are significant local variations between them. Our 
contention is that any scheme for representing shape must be 
able to perform, over and above interclass classification, 
finer classification within a class itself. In other words, 
we are concerned with the sensitivity of the representation 
to finer details. The problem before us then, is to attempt 
to achieve as fine a distinction between related shapes 
as possible within the framework of whatever representation 
we might choose to use for analyzing shape. Thus, in another 
and related sense, we shall inquire into the ability of 
a representation to capture the gross and fine features of 
a shape, and the degree to which they are decoupled from 


each other. 
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We shall, in this thesis, restrict ourselves to the 
domain of hand-written characters - as people would naturally 
write them. There are two excellent reasons for such a 
choice. For one, there exists a huge data base of hand-written 
characters collected by Karnick[4] and hence the time- 
consuming and tedious task of obtaining input data can be 
forgone. Moreover, we can match the results obtained against 
the judgement of a human observer, which in the final analysis, 
is what counts, 

1,2 DEFINITION OF SHAPE: 

It might be well to clarify our ideas as to what we 
mean by shape, and put it in mathematical terms. Since we 
are concerned only with 2D-shapes, we shall restrict ourselves 

o 2 

to the Euclidean plane, R . The Euclidean plane R is a 
two-dimensional vector space over the real numbers. Trans- 
formations of R^ which preserve all ratios of distances are 
called similarities. In R^, a similarity is either a trans- 
lation, rotation about some origin, multiplication by a 
constant of all distances from the origin (uniform scaling), 
or one of these followed by a reflection. 

We shsll start with the notion of a boundary or outline. 
An outline or boundary will be defined as a smooth closed 
curve in R^ i,e. a map of a part of the real line into R 
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which is twice differentiable 'almost everywhere' except 

at a finite number of points. Any outline in is mapped 

by the similarity transformations onto a four parameter 

family of other outlines. We can translate it to any origin 

(which is two degrees of freedom), rotate it (one degree of 

freedom) or alter its scale (one d,f.), two outlines are 

equivalent if one can be mapped onto the other by a similarity 

transformation in R , The similarity transformations form 

a group - the so-called equiform group - and hence we may 

construct equivalence classes of outlines which are collections 

of all the outlines that can be mapped onto one another by 

2 

similarity transformations, A shape in R can now be defined 
as an equivalence class, under the equiform group, of outlines 

O 

in R , Intuitively, we can consider a shape to be an outline 
which does not contain any information about position, scale 
and orientation in the coordinate frame, A shape repre- , 
sentation or measurement is defined as a function on some 
domain into the real line which is the same for all elements 
of an equivalence class. The domain of definition of the 
function can be specified in various ways depending on the 
nature of the shape variation to be studied. It can be for 
example, a finite point set, or the real interval, say (0,1). 
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1.3 SOME criteria FOR REPRESENTATION OF SHAPE i 

The most useful property of any representation for 
shape is that it can make some types of information explicit 
focussing our attention on the essential information, allowing 
us to look at the shape from another perspective which gives 
more insight. A shape representation, as discussed in the 
previous section, must be the same for all the outlines which 
belong to the same equivalence class. 

Another requirement that might be demanded of a shape 
representation is that it should contain information about the 
shape at varying levels of detail. Moreover, it should be 
clear from the representation which features of the shape belong 
to the coarser levels of detail and which to the finer levels. 
In other words, the degree of similarity between two shapes 
must be reflected in their representation, yet, at the same time 
even subtle differences should be expressible. These two 
conflicting requirements can be met only if it is possible to 
decouple stable information that captures the more general and 
less varying properties within a class from inf ormation tha t is' 
sensitive to the finer variations within the class* 

Finally, the representation should be capable of 


being computed easily. 
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V/e shall, at the end, attempt an evaluation of our 
scheme for representing shape - the so-called Fourier 
Descriptors - based on the present work and on [ 4 ] . 

1.4 ALGORITHMS FOR SHAPE ANALYSIS: 

There are many possible criteria according to which 
we can classify the methodologies used to analyze shape. 
First, we can distinguish between the algorithms which 
examine only the boundary points and ignore the interior 
area bounded by the boundary, and those which make use of 
the information contained in the interior region of a 
boundary as well, e.g. medial axis transformation (MAT). 

The boundary tracers can, in turn be classified into those 
which proceed sequentially along the boundary points; and 
thoee which examine all or most pairs of boundary points at 
each stage. Fourier descriptors and parsing of the boundary 
according to a simple grammar come under th€^ first category 
while parsing of the'o.:> boundary according to a context 
sensitive grammar belongs to the latter class. 

Again, we can distinguish between scalar transform 
techniques, which -transform a picture into a set of scalar 
features e.g., Fourier coefficient and space domain 
techniques which transform one picture into another, e,g, 
MAT. An even more fundamental classification is between 
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information preserving and information non-preserving 
techniques depending on whether it is possible to reconstruct 
a reasonable approximation to the original picture from the 
shape descriptors or not. The reader is referred to [ 8] 
for a more detailed discussion on the various shape analysis 
algorithms mentioned above, 

1.5 PLAN OF THE THESIS; 

Chapter 2 gives a brief review of the different 
approaches that have been adopted in the literature to 
describe a shape by Fourier methods. The methodology that 
has been followed by us is then presented in detail. In 
Chapter 3 we describe the two methods we imployed for 
classifying closely related shapes. One is based on the 
concept of bending energy and the other on that of matched 
filtering. Chapter 4 will deal with the experiments that 
have been done and the results obtained. 



CHAPTER 2 


SHAPE ANALYSIS AND FOURIER REPRESENTATION 

2.1 INTRODUCTION: 

In Chapter 1, we mentioned some of the schemes for 
analysing shapes. We shall present here the scheme that we 
have used for describing shape - the so-called Fourier 
descriptors. The method of Fourier descriptors constitutes 
an analytical description of shape in which the principal 
idea is to convert the description of a geometrical shape to 
an alternative, equivalent description by means of a linear 
operator. The ’spatial' domain description is replaced by a 
'frequency' domain description, • 

2.2 FOURIER DESCRIPTORS: A REVIEW:. 

Consider a plane closed curve C in the x ^y plane. 

Such a curve can be Fourier analysed in different ways. The 
Fourier descriptors of the curve are then defined in terms 
of the coefficients of the Fourier series expansion. It should 
be noted that this is an analytical method of describing 
a curve, an information preserving transformation of the 
original curve into a set of scaler features. Moreover, one 
can get as exact a representation as one wishes by simply 
considering sufficient number of coefficients. 
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The curve C can be represented by a set of parametric 
equations in the parameter t 

x=f^(t) 

y = fy(t) (2.2.1) 

The function f^^ and f^, are periodic in the parameter t and 
hence can be expressed in a Fourier Series. If we take the 
parameter t to be the arc-length s along the curve, measured 
from some arbitrary reference point on the curve, then we 
have a special form of the description of a curve parameterized 
by its own arc- length; 

X x( s) 

y = y(s) (2.2.2) 

Each of the functions x(s), and y(s) can be expressed in a 
Fourier series; 



I 

n=l 


a cos2 
n 


..2n7cs 

‘’ir~ 


+ 


n 


sin 


2 n‘jts 

L 


Y (=) = '=0 ^ Vi" ^ 


(2.2.3) 


where L is the total length of the curve. 

The set a ,b ,c ,d n=0,l, of the coefficients is 

n* n' n n 

defined to be the Fourier descriptor of the curve. This is 
the approach that has been followed in this thesis. 
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Gxanlund [3 ], Richard and Hemami [10], Persoon and 
Fu [9 ] among others have followed an essentially similar 
method as above but conceived of t'he curve as a complex 
function 0(is) 

0(s) = x(s) +r'^ y(s) 

The Fourier descriptors, F^* are the coefficients of the 
complex Fourier series expansion of 0(s). 

CO 

0(s) = S F < e 

n=_oo 

The curvature K at each point on the curve forms another 
important descriptor of the shape of a plane curve. The 
so-called intrinsic equation of a curve expresses the curva- 
ture K as a function of the arc length s, 

K = K(s) . , 

The parametric description of a curve, as in Fq. (2,2.1) 
and implicit description of the curve as f(x,y) = 0, contain 
information about the position, and orientation of the curve 
relative to the coordinate axes as well as information about 
the shape of the curve. This is reflected directly in the 
frequency domain ,by the dependence of the Fourier coefficients 
on the position and orientation of the curve in a coordinate 
system. The intrinsic equation of the curve describes only 
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the pure shape of the curve without reference to any coordinate 
axes. Since, moreover, the intrinsic equation of a closed 
curve is periodic in s i.e., K(s+L) = K(s), we can expand it 
in a Fourier series. Unfortunately however, for curves which 
have corners, the curvature is infinite at the corner points 
and the method fails. The tangent angle 0 at any point on the 
curve is related to the curvature K through the equation: 

^ ~ ds 

If we define the cumulative angular function ©(s) as the 
difference between the tangent angle at the point s and that 
at the starting point s= 0 , then we have 

0 (s) = 0 (s) - 0 ( 0 ) = /\.ds 

0 

* 

Defining © (t) as 

e^(t) = 0(|~) 4- t 0 < t < 271 

we get a function which is invariant under translation, 
rotation and scale change. Zahn and Ruskies [l2 ], Rennet and 
MacDonald [ 2], among others, expand ©*( t) in a complex 
Fourier series. 

©*(t) = 2 

n=~°° 
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and the set n = .1,2, , . forms the FouriE 

descriptor of the curve. 

It is interesting to compare the different methods 
of obtaining the Fourier descriptors. The function 0*(t) 
in the tangent angle method is almost always discontinuous 
and consequently the Fourier coefficients can decrease no 
faster than c/n, where c is a constant independent of n, as 
n ®°, On the other hand, the parametric functions x(s) and 
y(s) are always (piece-wise) continuous functions of s; 
hence the coefficients in their Fourier expansion decrease at 

O 

least as rapidly as c/n as n-^^. Another important difference 
is that a subset of the descriptors F^ usually does not give 
a closed curve, whereas a subset of F^* does. One more 
advantage of using the coordinates x(s), y(s) instead of the 
tangent-angle representation of the boundary curve is that the 
parametric representation is less sensitive to the noise 
inherent in a fuzzy boundary. The tangent angle is related 
to the derivatives of x(s) and y( s) . A direct consequence 
of this is that small variations in the coordinate values of 
the boundary points can result in large variations in the 
tangent vector. 
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The Fast Fourier transform or direct integration has 
oeen used in the past to calculate Fourier descriptors. 
However, if the contour is represented by a piece-wise 
linear approximation, a more direct method of calculating 
the Fourier coefficients can be employed, as has been shown in 
Kuhl and Giardi’na U 5 ] who also have described a normalization 
procedure, which we will borrow, based on the elliptic 
properties of the Fourier coefficients. As mentioned before, 
the Fourier coefficients are dependent on the choice of the 
starting point of the traversal round the contour, and on 
translation, rotation and uniform scaling of the contour. 

In many applications in pattern recognition, the 
position size and rotation of an object are not important. If 
we describe the. object by its boundary, the starting point 
is irrelevant also. Hence, self-consistent normalization 
procedures based only onfhe intrinsic properties of the curve 
must be used, 

Granlund [ 3] considered products of Fourier series 
coefficients which are invariant to position size, orientation 
of the contour and starting point choices * The result is 
an increase in data dimensionality from n to n/2 without 
any change in total information, persoon and Fu [9 ] deal 
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with unnormalized coefficients and reduce the problem of 
shape matching to one of finding an optimum size, orientation 
and starting point match between the sample FD (Fourier 
descriptors) and each reference FD in the training set. 

Although this method retains all the shape information of 
the original contour, and an optimum match for size, orientation, 
and starting point can always be found, it is computationally 
expensive, Richard and Hemami [10] follow a similar matching 
criterion in their work on identification of airplanes. 

The invariant elliptic properties of the Fourier 
coefficients can be exploited for a convenient and intuitively 
pleasing scheme of normalization, so as to get pure shape 
information, discarding irrelevant information introduced by 
the rotation, translation, size and starting point factors. 

2.3 THE J\£TH0D: 

2.3.1 Fourier Descriptors of a Fiece-wise Linear Contour: 

Consider a plane closed vurve C in the x-y plane which 
is described in terms of the parametric equations ■ 

X = x(s) 

y = y(s) 

where s is its own are- length measured from some arbitrary^^ ^ ^^^^ ^ ^ 
reference point on the curve. The parameter s varies f rom ^ ^ ^ 
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0 to L where L is the total length of the curve. The x and 
y projections are periodic functions in the parameter s, 

X (s+L) = x(s) 
y (s+L) = y(s) 

We can express x(s) and y(s) in a Fourier series: 


x(s) = A +• E a cos + 2 b sin 

° n=l ^ n=l ^ 


y(s) = C + S c_ cos + E d^ sin 

° n=l ^ ^ n=l ^ 


The coefficients a ,b »c and are given by 

n n n n 

, L 1 L 

Ao = L / x(s)ds, c^ = Y zf y(s)ds 


(2,3.1. 1) 


o 


'0 “ L 


= a / x(s)cos ap ds, = a / ^(3) cos ap 


ds 


= 1 f x(s) sin ds, = f / y(s)sin 


2n7ts 


ds 


(2.3. 1.2) 


Now suppose that we use a piece-wise linear represen- 
tation for x(s) and y(s) i.e. we consider C to be described by 
a set of discrete points (Xjj^*y^) i= 0,1, * . .m , Then, the^^^ 
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derivatives and ^ consist of a sequence of piece— wise 

constant derivatives and Ay / as in the interval 

P P P P 

®p_l < s < Sp as p varies from 0 to m. We note that the 
cix 

derivative ~ is itself periodic and hence can be expanded in 
a Fourier series ; 


dx 


2n'n;s 


. _ 2mts 




where 


^ _ 2 fT dx — 2mts 

“n - L J Hi 


_ 2 p . 2n'n:s 

- L J di “r“ 


Thus, 


a = 
n 


2/42 cos 2 £Ss ds 


L ds 


2 « ^ 


— y 
L ^ As„ 
p=o p 


/ cos 2p ds 
®p-l 


a_ ” 4^(sin - sin 3l!S|E=i 

™ p=o -"'p' I- L 


n 


2 n dx ^ * 2nTcs . 
r / ^ Sin sds 

L J ds L 


1 I -lie /"^ 

L ^ <4s ■i. 

P=° P ®n-l 


2n'n:s 

sin —j”” ds 
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1 m ^ X 

— — y Hi 

nit ^ s. 


( cos 


2n7is. 


2mts^ 1 

— ) 


H V 

Now, we can obtain by directly differentiating the 
Fourier expansion of x terra by term. 


y _ 2mt ^n-Fus 2mc ^ 2mrs 

1 - L \ "L “ ^n "1“ 


Equating coefficients obtained in the two expressions for 


dx ^ . 

^ we get: 


T ni AX 2n-jts^ 2mts , 

= -jV- [co= ---2 - cos ^E=i ] 

2n -n; p=o p ^ 


, m X ^iims ^mis^ , 

2-r ^ -Kf [cm -T-a - Sin 

2n^ic"' p=o ^ ®p ^ ^ 


2mis 


2mis 


Similarly for the y-projection we get 


r m Ay 2nits 2niis . 

c^ = [cos ^ coS ] 

" 2nV p=o -^Cp I- ^ 


j m Ay„ 2nits 

2n Tc p=o p 


2^13^ 


2,3,2 Elliptic Properties of the Fourier Coefficients: 

The Fourier series representation for the x and y 
projections can be written as: 
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N 


x(s) 

II 

> 

z 

O 

n=l 



N 

y(s) 

II 

o 

> 

+ 

Z 



n=l 


o < s < L 


where, 


Xn(s) 

Y„(s) 



COS 

cos 


2n7cs 

L 

2mis 

L 


* '•n 


sin 

sin 


2mis 

L “ 

2n2S 
L “ 


(2. 3.2.1) 


( 2 . 3 . 2 . 2 ) 


For a given n, .the points form an elliptic locus 

as s varies from 0 to L. After a little algebra with 
equation (2, 3, 2, 2) we get 


(d ^ 4 -c ^)X ^+(a ^)Y ^-2(a c +b d )X Y 

^ n n ^ n n n ^ n ^ n n n n^ n n 

(ad - b c )^ 

' n n n n^ 


= 1 


which we recognise to be the equation of an ellipse. The 
fact that for each n the points describes an ellipse 

allows us to interpret in an interesting way the Fourier 
series representation of a contour, V\/e can conceive of the 
original contour as being traced out by the phasor addition 
of rotating phasor s, one for each harmonic. Each phasor 
sweeps out an ellipse, and the phasor associated with the 
nth harmonic rotates n times, as fast as the fundamental. 

More important from the point of view of normalization 
is the fact that the elliptic loci independent 
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of the starting-point of the trace round the contour. To 
see this let 


Then we have 




s'+s^) 


from equations 

2n'n;(.s* +S q) 2n-n;(s'+s ) 

a cos 4 . sin -f 

n L n L 


Yn(s’+So) = cos 


2nit(s'+s^) 


■r— - — + sin y 
L n L 


2mi(s’+s^) 


This can be rewritten as 


(2. 3.2. 3) 


= V V 




where we have put 


X^‘(s*) = X^(s‘4.s^) 


(2. 3.2.4) 


Y « (s' ) = Y„(s'4*s^) 
n ' ' n^ o 


(2. 3. 2.5) 


The coefficients a^' , b^' , c^' and d^' are the 
coefficients of the expansion for the starting point s'=0 
i^e., a starting point shifted s^ units along the curve from 
the original starting point s=0. The new set of coefficients 

a ',b ',c*, d„' can be obtained from the old set 
n ^ n ' n 

a ,b ,c ,d through the transformation T . 
n’ n' n* n 
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a ' c • 
n n 


COS 

2nTCS 

0 

L 

2n-ji: ^ 
Sin s 

L 0 


^n 


b„* d • 
n n 


-sin 

2n-n:s 

0 

L 

^^^2n'n; 

cos r— 

L 0 



'‘n. 



- 






(2*3. 2. 6) 

Again, eliminating sines and cosines from equations (2. 3. 2. 4) 
and using equations (2. 3.2, 5), (2, 3. 2. 6) we get 


(d„^-l-c^^)x‘^ +(a^^+b ^)y‘^ - 2X »Y ’ (a c +b d ) 
n n n ^ n n ^ n n n ' n n n 


= 1 


(ad - b c )‘ 
'■ n n n n'^ 


which gives the same locus as before. Hence the same elliptic 
loci are obtained for different starting points. It is this 
invariance that we will exploit for the starting-point 
normalization. 

2,3,3 Rotation of the Contour: 

Any rotation of the contour, with respect to the 
coordinate axes is reflected in the Fourier coefficients. To 
see the effect of this rotation on the coefficients, assume 
that the X,Y axes are rotated through a counter clockwise 
angle 0 into a new set of U,V axes. The new coordinates 
will then be given by the familiar transformation: 
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u 


COS0 

sin0 

. V 


-sin0 

COS0 _ 


The and projections on 
in matrix form as 


X 

„ Y_. 

the X and Y axes may be written 




s 

n 

b 

n 


cos 2"^® 






sin 2™= 


The projections U^, on the U,V axes will then be 


Un 

' COS0 

sin0 


' V 


''n 

~sin0 

COS0 


.""n. 



COS0 

sin0 

a b ” 

n n 

t-cos 


-sin0 

COS0 

c d 

L ^ 

_sin 2235= 


~a ’ ' 
n 

b ’ ’ 1 
n 


COS 

2^15 I 

L 

c ' ' 

_ n 

c 


sin 

2n‘jts 

L J 


where the axially rotated set of coefficients * ^n' ’ » 


Cn' ' , ’ are given by 
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a ‘ ' 
n 

V 


COS0 

sin0 


a 

n 


c • ' 
n 

V 


-sin0 

COS0 


„'=n 



Finally, if we shift the starting point by s^ units 
and also rotate the X,Y axes by 0 degrees in the counter- 
clockwise direction, then the combined effect on the Fourier 
coefficients can be expressed concisely as 


a > ’ 
n 

b • • 
n 


COS0 

sin0 


"a 

n 

b 1 
n 


2nTts 

cos ° 

2mts -1 

- i „ 0 

-sin 

c ' » 

n 

b ' ' 
n 


,j“Sin0 

COS0 

1 

c 

L n 

d 

n . 


2nTts^ 

[sin 

2n%s^ . 
0 

cos ^ 


2.3.4 Uniform Scaling: 

Consider the expression for the coefficient a^ 

= I A(s) cos ds . 

o 

If we put = 9 we get, 

i 2n; 

a=“-/x(9)cosn9d0 

nut 

0 

If we blow up or ■. ..x the original contour by a factor K, 

we see immsdiately . from the last equation that the coefficient 
is multiplied by the factor K, Hence, scaling the contour by 
a factor K scales all the coefficients by the same factox K. 
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2.3.5 Normalization of the Fourier Descriptors: 

Vve have seen in the last section that the Fourier 
coefficients contain information about the starting-point, 
spatial rotation, size and translation of the contour 
relative to a coordinate system, as well as information 
about the intrinsic shape of the curve. Hence normalization 
must be performed to extract the pure shape information. The 
invariant elliptic properties of the Fourier coefficients 
form the basis for a convenient method of normalization. 

The coefficients and contain the d.c, infor- 
mat ion or the translation of the contour in the reference 
frame. By simply ignoring these terms we can achieve 
independence with respect to translation. The other factors 
pose a slightly tricky problem. The normalization is done 
in two stages. To start with, the first harmonic phasor is 
rotated until it is aligned with a semi-major axis of its 
locus. This normalizes the starting point. Then the X,Y 
coordinate axes are swung into a new set of U,V axes 
defined by the major and minor axes of the ellipse such 
that the positive X axis is coincident with the semi -major 
axis located above. This achieves normalization with 
respect to rotation of the contour, ■ 
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Consider the first harmonic locus 


where Q. 


^1 ~ + bj^sin© 

Yf = c^cose + dj^sin© 
2tis 


2n'n:s 

The starting-point shift s^, or angular shift — ° can 
be calculated by finding the value of © for which the 
magnitude of the first harmonic phasor E = is a 

maximum. Differentiating E with respect to © and equating 
it to zero^ gives us 


e. 



^^^ 1^1 ^ ^ 1 ^ 1 ^ 

37^2 vTT"2 
a i"^c 1 — bj^— dj^ 


That this gives a maximum is checked by evaluating ^ 

d© 

at 9 = 0Q and noting that a negative value is obtained. 

The Fourier coefficients, a^‘ , , d^*' for 

the displaced starting point are given by 





cos©^ 

sin©„ 


o 

0 


_-sine^ 

cos©^ 


v/ 




To determine the spatial rotation 0 of the contour we write 
the first-harmonic locus of the point 1 
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h' 

(s') = 

cos 


sin 

'^1 

(s' ) = 

cos^pl 

+ 

sin 


Since when s’=0, the first-harmonic phasor is aligned with 
the semi-major axis, the angle 0 is given by 

0 = tan- [ 3^^^] 

= tan”^ [ ] 

^1' 

The G axially rotated set of coefficients will then 
be given by: 


a f ’ 
n 



COS0 

sin0 


a ’ 
n 

b »" 
n 

c » ' 
^n c 

d ' ' 

n 


-sin0 

COS0 


c ' 
n 

-4 


If we divide each of the coefficients by the length of the 
semi-major axis V” ( o)^+Yj_' (o)^ = f (a^* )^+( then 

the coefficients are made independent of size as well. 

2.3.6 Some Remarks; 

The Fourier series representation for shape requires a 
closed contour. Open outlines have been effectively made 
closed by retracing the outlines from the end back to the 
beginning. ?i/e are currently working with 12 harmonics as this 
gives a fairly close (2 pixels) approximation to the original. 



CHAPTER 3 

CLASSIFICATION OF SHAPES USING BENDING-ENERGY AND 

MATCHED-FILTERING 

3.1 IINTRODUCTION: 

An intuitive way of classifying shapes is on the basis 
of how ' Wiggly* they are. Translating this intuition into 
mathematical terms, we can employ curvature or the rate of 
change of the tangent with respect to the arc length, as a 
descriptor of shape. Curvature by itself and measures based 
upon curvature have played an important role in shape repre- 
sentation and pattern recognition. Indeed, it will not be 
an exaggeration to say that curvature, either directly or 
indirectly has played a prominent part in most shape analysis 
schemes. If we are to believe some of the theories of vision, 
curvature plays an important role in human perception of 
shape [ ]. Unfortunately^ curvature is too sensitive to 

noise inherent in a fuzzy boundary to be of much use as a 
shape descriptor when used by itself. On the other hand, if 
we average out the curvature over a length of the boundary, the 
average might serve as an useful measure of how bent or 
curved the boundary is Interestingly enough, the stored free 
f' energy in an elastic bar under deformation Is proportional 



28 


2 

to /K dt where K is the curvature and t the length along the 
bar. Such measures have been used to classify biological 
shapes [ ll] , A shape, seen from this point of view is a 
deformed rod. The bending energy of such a shape can be 
related^ to the Fourier coefficients in a simple manner as 
will be shown in the succeeding section. 

Another idea for classifying shapes, based upon energy, 
that we have pursued is that of optimum filtering or the 
so-called matched-filtering, Vve can consider eachs shape as a 
signal and the problem of shape classifications, becomes' one of 
determining how closely the signals clutter together in the 
signal space, 

3.2 CLASSIFIC AXIOM OF SHAPE USING BENDING ENERGY: 

The shape that an elastic bar assumes, when its ends 
are constrained, is determined by the condition that the free 
energy associated with the deformation shall be a minimum. 

If it is assumed that the cross-section is small in comparison 
with the length of the bar, and further that there is no 
appreciable torsion of the bar, then the bending energy is 
proportional to /K^dt, Consider a rod of negligible cross- 
section deformed as shown in Fig. 3,2.1. Then the free energy 
stored in the rod at a point t on it is given by [6 ] 
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e = I Y(t) Iy(t) K^(t) 

where , 

Y(t) is the Young’s modulus at the point t, 

I^(t) is the moment of inertia of cross-sectional 
area about the y axis, and 
K(t) is the curvature at point t. 

The total energy over the rod of length T is 

£ = I / Y(t) Iy(t) K^dt 

If the cross-sectional dimensions are negligible compared to 
the length of the rod, then Iy(t) is a constant, say I; if we 
assume further that the Young's modulus is also a constant Y, 
then we have 

T 

E = ^ YI / K^(t)dt 
o 

If we put =E/(r/ 2 ^YI) we get 

E_ = f K^dt (3,2,1) 

0 

Now consider a plane closed curve (Fig. 3.2.2) 
param.eterized by its own arc length t, Vve shall adopt the 
vector notation as it is more convenient. Vectors will be 
represented by a symbol with an arrow above, and unit vector 
will be identified by a cap above the symbol. Referring to 
Fig. 3.2.2, the position vector R is given by 





X 


• / X - / y 

where a and ^ are unit vectors along the x and y axis 
X y 

respectively. The unit tangent vector T will be then given 
by 


T = 


dR 

dt 


dR . dx ^ ^ 

■ dt ®x + dt ^y 


The rate of change of the tangent vector with respect 
to the arc length is 


dt 


d X A , d^ A 


'71 “"x 


dt 


dt 


o a 
±.2 y 


But we know that, 


dl 

dt 


dT 

d^ 


dt 


= n K(t) 


where. 


n 


K(t) 


dT 

= ^ 



is the unit normal to the curve 
is the curvature at point t 


Hence, 

1 = K^(t) 

_ ( d X ^2 ^ v^2 (3,2,2) 

dt-^ dt^ 

Substituting in equation (3*2,1) we get. 
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= / [(^ 2 )^ +■ ]dt (3,2,3) 

^0 dt"^ 

Now, consider the Fourier series expansion for x(t) and 
y(t). 


x(t) = + Z cosnw^t + sin nw^t 

n=l 


y(t) 


00 


^o ^ 


Z 

n=l 


where 

o 



c cosnoa^t + d sin nw t 
non o 


(3.2.4) 


Parseval's relation states that for any signal f(t) 
represented by a Fourier series 

00 

f(t) = A + Z a cosnw^t -t- sin nw^t 
° n=l 


the following relation is true; 


T ^ o o 

f^(t)dt = + -I S (a^ + ^n 

o n=l 

From eqns. (3.2.3) and (3.2.4) and using Parseval's relation 
we get, 


T, I 2 (nco )'^ ‘^''n +^n ^ 

n=l 


(3.2.5) 


If we define an average bending energy per unit length 
then we get 
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E~ 

E — — i 

N ~ T • 

OO 

= 2 ( 3 . 2 . 6 ) 

n=i 

A little clarification with regard to the last equation 
is in order. If we are to leave the coefficients a^,b^,c^ and 
d^ unnormalized, then the bending energy will not be the 
same for an outline and an uniformly scaled version of it. If, 
for example, we blow up the outline by a factor m, the magnified 
outline will have a bending energy that is lesser by a factor 
m^. Although this observation is what we should intuitively a 
expect, yet we need a shape factor that is identical for all 
boundaries in the same equivalence class. On the other hand, 
if we choose to use normalized coefficients, we get a shape 
factor that is the same for all related, translated, and 
scaled versions of the original outline. The normalized 
coefficients, thus, fit in very neatly in this scheme. 


Equation (3.2.6) reveals some very interesting features. 

The bending energy is proportional to a weighted sum of the 
squares of the Fourier coefficients. The presence of the 
factor n"^ implies that as n increases, the weight given to the 
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harmonic coefficients increases very rapidly. The higher 

harmonics are given far more importance than the lower 

harmonics. For example, the relative importance .ascribed 

to the first harmonic and the tenth harmonic will be in the 
4 

ratio IslO-, Now, we intuitively expect that the lower 
harmonics determine the gross features of the shape while the 
higher harmonics control the finer details. Pattern recog- 
nition work done by Karnick [4], Fu [9] and others provide 
experimental evidence to support such an hypothesis. Keeping 
this fact in mind, then, our scheme of weighted sum with the 
weights increasing exponentially with n seems well suited for 
distinguishing between shapes whose general pattern is 
unvarying but which exhibit subtle differences between each 
other. It will be seen in the next chapter that such indeed 
is the case. In such a classification, all shapes with 
identical stored energy in their contours will belong to the 
same equivalence class. 


It is possible to make still finer classifications of 
shapes by considering the energy stored over small lengths 
of the boundary, instead of the total energy over the whole 
permiter. Indeed, the incremental aergy stored in the. 
boundary, between tj^ and t^ will be 


A E 



K^(t)dt 


(3.2.7)! 



36 


From equations (3.2.2), (3.2.4) and (3.2.7), we get, 

AE = / ^ [(^)2 + (^)^]dt 

tj^ dt dt 

^2 N 2 2 

= / ( 2 (nw ) [a^cosnto t + b^sin na)^t])'‘'^dt 

+.*' 'l O ^-n o n O'" 

tj_ n=i 

N 2 2 

+ ( £ (nu) ) [c^ cosnw t-i^d sinnco t])^dt 

o n o n o"^ 

n=l 

(3.2.8) 

The RHS of equation (3.2.8) can be integrated using 
numerical methods. 

To illustrate the usefulness of such an r'.nv incremental 
measure suppose that we are required to distinguish between a 
shape and a reflected image of the shape (see Fig, 3,2.3). The 
total energy in the shape will be the same for the shape and 
its mirror image. However if both the contours are traversed 
in the same direction, assuming that we start the traversals 
at the points marked, then, it is obvious that the energies 
stored in the lengths ab and a^b' will be different, 

,3.3 CLASSIFICATION OF SHAPE USING MATCHED FILTER: 

We can consider the issue of recognizing a shape as 
that of detecting a signal in the presence of additive noise. 
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Fig, 3,2, 3t A shape and its laixror image 
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In detection, one is concerned not so much with reproducing 
the original signal as with being able to establish the 
presence or absence of the signal. Typically, the form of 
the signal is known in advance and one is required to decide 
whether the signal component is present in the receive'd wave, 
corrupted by additive noise. It turns out that when the 
power spectral density of the noise is constant i.e, when the 
noise is white, the optimum solution to the problem of 
detection is a matched filter. The detection system in such 
a case is called a matched filter because its characteri- 
zation is matched to that of the signal component, in the 
received complex of signal plus noise, A matched filter is 
optimum in the sense that it maximizes the signal to noise 
ratio at the filter output; this ratio depends only upon the 
ratio of the signal energy to the spectral density of the 
white noise at the filter input. 

The Fourier series representation for x(t) and y(t) 
suggest that the signals x(t) and y(t) can be considered to o. 
be vectors in the infinite dimensional vectox space in which 
the set cosnco^t, sinnw^t; n = 1,2, . . . form an orthogonal 
basis set. We know, however, that in practice all signals: 
must have a finite energy and hence thecoef ficients a^ and| b^ 
must tend to zero for large values of n. Hence, although the 
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dimensionality of the signal space is infinite we can, for 
all practical purposes consider an approximate finite 
dimensional space with dimension N, Any signal, then, can be 
considered to be a point in this N~dimensional space. Viewed 
in this perspective, optimum filtering turns out essentially 
to be a scheme wherein, given a set of m signal points 

s^(t) i = 1, ...m in the signal space, the unknown signal 
is declared to be identical to that signal for which the 
distance from the unknown signal is a minimum [ 7]. 


The received signal s^ is declared to be that signal 
Sj^, for which the decision function 



I 

n=l 




oo 


) “ 2 
n=l 



a 


rn 



(3.3.1) 


is a minimum 


It should be noted that when i=r, D=o. 



CHAPTER 4 

EXPERIMENTS^ RESULTS AND DISCUSSION 

4.1 INTRODUCTIONt 

As we stated in the beginning, one of our main aims 
was to investigate the degree to v\tiich the Fourier series 
representation fox shape can achieve decoupling, in the 
representation, between the gross and fine features of a 
shape and we stated in the previous chapter oux belief that 
the lower harmonics dictate the gross features of the shape 
and the higher harmonics the finer features. In Chapter 3, 
we also argued out how, based on the above assumption, the 
bending energy associated with a shape constitutes a good 
basis for sorting out variations within a class, l/ve also 
pointed out how classifying shapes within a class may be 
considered within the framework of matched filtering of 
signals, 

4.2 THE GLOBAL VERSUS LOCAL NATURE OF THE FOURIER REPRE^NTATION 

As seen from a purely mathematical point of view, the 
Fourier coefficients have the property that of all trigono- 
metric polynomials of order n, that polynomial which has as 
its coefficients the Fourier coefficients of the function f(x)p 
has the least mean-square deviation from that function f(x). 



41 


The Fourier series of a function f(x) converges to the function 
f(x) at all points where f(x) is continuous and to the mean 
value of the jump where there is a jump discontinuity. Since 
the coefficients of the Fourier series are obtained by inte- 
gration, any change in any part of the curve will affect all 
the coefficients. While, it is true that a drastic change 
such as a step discontinuity introduced over a small portion 
of the curve, will adversely affect the convergence of the series 
everywhere (the series will now converge more slowly), any 
small change in the function values in a narrow neighbourhood 
which does not destroy the differentiability of the function 
will affect the higher harmonics more than the lower harmonics. 

To what extent our assertion that the lower harmonics 
control the course of the function in a global sense while 
the higher harmonics are responsible for local variations can 
be seen from figures 4,2.1 and 4,2.2. Strokes, lobes, curls 
etc. are regions of high curvature and contribute to the 
presence of significant higher harmonics in the Fourier 
representation. 

Part (a) of Fig, 4.2.1 shows a shape generated with 
all the twelve coefficients being taken into account. Part (b' 
shows the figure obtained by considering only the .• first 


four harmonics. 
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(a) (b) 


Fig, 4.2.1: Illustrating the role of lower and higher 
harmonics, 

a) Shape re-generated with all the 12 harmonics 

b) Shape regenerated by considering only the 
first four harmonics. 
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Fig. 4.2.2: Illustrating the role of lower and higher 
harmonics. 


a) Shape regenerated with all the 12 harmonics 

b) 8-harmonic approxiiretion 
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It is evident from the above examples that our 
contention that there is a decoupling, fuzzy though it be, 
between the gross and fine features of a shape, in terms of 
the lower and higher harmonics, is quite admissible. 

Moreover, Karnick [4 ] has found that reliable Interclass 
separation of hand-written characters can be achieved by 
considering only the first five harmonics. This corroborates 
our claim that the more stable and less varying properties 
of a shape are captured by a few lower harmonics. 

4.3 CLASSIFICATION BASED ON BENDING-ENERGYs 

In Section 3.2, we savj how a boundary or outline of an 
object can be considered to be an elastic bar undergoing 
deformation. We also saw how the free energy stored in the 
shape can be related to the Fourier coefficients via equation 
(3,2.6). This point of view leads to the use of the bending 
energy associated with a shape to generate a set of equivalence 
classes, all shapes within a class being equivalent in terms 
of the stored energy in them. 

Fig. 4,3.1 shows three shapes which are closely allied 
to each other and the bending energy associated with them. 

It should be noted that the classification is in harmony with 
the way an human observer would classify them, dc ' based on 
his intuitive ideas of how curved ’or’ 'wiggly' a shape is. 
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B.E. = 288,805 


B.E. = 530.897 


B.E. = 766.118 


Fig. 4.3.1: Three shapes and theii bending 
energies. 
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Some more shapes and their bending energies are 
given in Fig* 4,3*2* 

Fig* 4,3.3 shows fifty one variants of the character 
’C’ arranged in the descending order based upon bending 
energy* The values of the bending energies are listed in 
Table 4.3,1. 

We see from Fig, 4,3*3 how the smoother variants — those 
devoid of any strokes, lobes, curls etc, cluster together 
at the lower end of the energy spectrum and how the more 
Wiggly ones occupy the higher end of the spectrum. It should 
be pointed out further that shapes which are quite close to 
each other are also classified at the same level. For 
example, the variants C<22>, C<26> and C<27> are quite close 
to each other and this is directly reflected in terms of 
their bending energies. It is also a vindication of the: fact 
that the Fourier series representation is stable - the degree 
of similarity between two shapes is reflected in the repre- 
sentation as well, which is sufficiently insensitive to minor 

variations. 

The results of similar classifications performed on 
the set of 51 variants of the character W and S are given in 
Tables 4.3.2 and 4.3.3 respectively. The reader is referred 
to the Appendix for a list of the figures. 
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B.E, = 83,838546 


B.E. = 70,259707 


B.E. = 465.14971 


Fig, 4.3,2: Bending energies of three variants of B 





C<30> C<i4> C<15> C<36> 



C<10> 



CO> C<35> 



•C<28> 




C<42> 



C<33> 


Fig. 4«3.3t 51 ysriants o# ‘c* azi'anged in D^«e«nding oid«r 
of Bonding Energy* 


C<21> 


C<3^> 


C<7> 


C<32> 



c<:i7> 



c<o> 


C<37> 



C<49> 




C<16> C<^1> 


A 

A 

C<8> 


/ 


C<23> 



C<24> 


C<26> 



C<22> 


C<27> 


4#3*3t c^ntiiiiied# 





C<29> 



C<i2> C<:25> 





Cc6> C<47> 


r 

/ 


C<18> 



C<48> 
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C<38> 
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Table 4,3,1 


The bending, energies of 51 variants of ' C 

descending order 


arranged in 


Variant 


Bending Energy 


1 

43 

45 

2 

30 

14 

15 
36 
10 

9 

35 

31 
28 
20 
42 

33 
21 

34 
7 

32 
17 

0 

37 

49 

16 

41 


160.75951 

88.822914 

86.054553 

83,840840 

80.947044 

69.530754 

69.527396 

66.638730 

61.981963 

59.903209 

59.690878 

54.106864 

52*050690 

50.768181 

49.461648 

45,966666 

45,281404 

43.955914 

41.664268 

41.302578 

40.946887 

37.394480 

34.373846 

32.881071 

32.221905 

31*470717 
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V aEian t 

8 

23 

24 
26 
22 
27 
19 
44 

5 
13 
40 
50 
11 
39 
29 
12 

25 
4 

6 

47 
18 

48 
38 
46 

3 


Bending Energy 
29.725692 
29.586152 
29.441790 
28.719499 
28.086087 
27.909746 
26.613102 
25.322420 
25.186769 
23.717210 
21.864629 
21.689877 
21.093493 
19.772351 
19.004930 
18.531553 
18.424521 
17.905297 
17.671253 
17.572102 
17.162462 
16.792473 
15.712571 
14.330127 
14.270069 



Table 4.3.2: The bending energies of 51 variants 
of ’ i/l/’ arranged in descending order 


Varian t 

45 
0 

12 

49 
15 
48 

46 
3 

40 
38 

36 

50 
9 

11 

7 

18 

21 

28 

37 
27 
14 
25 
23 
10 

41 
34 

47 


Bending energy 


805.92713 

571.74282 

383,86132 

382.55382 

367.47211 

294.20727 

286.87660 

229.65457 

207.92111 

205.52481 

199.04057 

183.43976 

181.18240 

176.17072 

169.44450 

168.02066 

163.79804 

162.33990 

160.06561 

158.76482 

155.52497 

152.86746 

152.46783 

152.15717 

149.29733 

144.04170 

143,48608 
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Variant 

Bendinq enerav 

44 

124.36388 

13 

117.90828 

35 

116.69977 

33 

107.86908 

29 

104.77307 

19 

102.62156 

39 

98.361093 

30 

94.112065 

32 

90.154442 

1 

89.851562 

42 

88.846672 

17 

86.225944 

20 

82.617205 

26 

79.992243 

16 

78,948387 

2 

78.250048 

5 

72.891058 

22 

69.170465 

43 

66.820249 

6 

64.118934 

31 

59.640653 

8 

58.458732 

4 

56.992058 

24 

44.251450 
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Table 4.3,3 

The bending energies of 51 variants of 'S' arranged in 

descending order 


Variant 

42 

13 

39 

40 
16 

37 
15 

38 
17 
47 
46 

1 

2 

49 

44 

43 

41 
0 
7 

32 

30 

45 

14 

50 
3 

29 

5 


Bending energy 

660.73912 

159.05835 

142.65612 

140.84575 

129.67440 

122.98946 

122.71067 

113.83172 

86.963529 

85.963388 

81.467459 

31.042984 

80.037961 

69.380047 

68.725529 

60.638382 

59.256546 

58.791855 

54.907763 

53.826262 

53.520165 

51.798248 

51.957600 

44.494440 

44.239019 

42.548772 

39,555336 



Variant 


36 

31 

35 

4 

10 

27 

19 
6 

26 

28 
9 

18 

22 

11 

21 

12 

25 

48 

8 

33 

34 

20 
24 
23 


Bending enemy 

39.045345 

38.806584 

38.039027 

36.563881 

35,033606 

34.559309 

34,180757 

33.972211 

31,898164 

30,677968 

30.003176 

29.004076 

28.521631 

28.394356 

23.058547 

22.806727 

21.697493 

21.649204 

21.421442 

19.400690 

18.255912 

15.821547 

14.628535 

14.524376 
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It is possible to make finer classifications if we 
employ the strategy of computing energies stored in smaller 
lengths of the boundary, as discussed in section 3,2, instead 
of the total energy over the whole length. To explain the 
philosophy of the method, we pose the problem of distinguishing 
between a shape and a mirror image of it. Obviously, a shape 
and a mirror image will have the same total bending energy. 

Fig, 4,3,4 shows a shape and its mirror image. The 
energies stored in the corresponding parts of the two shapes, 
calculated using eq, (3,2,7) is given in Table 4.3,4, It is 
evident from Table 4,3,4 that although the total energy in 
the two shapes are the same, yet a distinction can still be 
achieved on the basis of energies associated with corresponding 
parts, , 

The drastic reduction in the dimensionality of the 
data from n to 1 makes this method too insensitive to be 
suitable for obtaining reliable classifications in the larger 
space containing all hand-written characters. However, if we 
restrict our attention to the subspace containing the variants 
of a particular character, then the bending energy measure 
provides an intuitively appealing basis for discriminating 
between closely related shapes* 
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Fig. 4.3. 

Table 4,3 

Original : 

Length of 
curve 

ab 

be 



4: A shape and its Mirror Image. 


.4: Bending energies of a 


shape 


the 


Bending Energy- 


shape and its mirror 

Mirror Image 
Length of the curve 


288.80488 4 2- 
383.05923 


image 

Bending energy 

383.05923 
(288.80488') T 2- 

£527.46167 


Total 


527.46167 


Total: 
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4.4 CLASSIFICATION BASED ON MATCHED-FILTERING: 

In this approach, we choose a prototype and decide 
whether a given shape is an instance of the prototype by 
testing if the decision, function D given in equation (3.3.1) 
is zero or not. In order to classify shapes we assume a 
certain threshold, fixed at a certain percentage of the energy 
of the prototype, and put those shapes for which the decision 

function D is less than the threshold in the same equivalence 
class. 

The first prototype we chose was C<2>. All the 
51 variants of the character C were passed through the filter 
matched to the prototype C<2>. Table 4.4.1 shows the classes 
obtained when the threshold was fixed at lQ/4 of the energy 
of the prototype and-at 20yi» 

The next prototype we chose was C<3>. The results 
obtained in this case are shown in Table 4.4.2, ' 

- * 

The prototypes were chosen based on the results of our 
experiment 'using the bending-energy concept. A look at table 
4.3.1 shows that these two prototypes are widely separated in 
their energy spectrum - C<3> lies at the bottom of the range 
while C<[2> lies at the top. This observation leads us to expect 
that the sets of the equivalence classes obtained in the two 
cases should be ■Jisjoint, Vife find that such indeed, is the , 
case from Tables 4.4.1 and 4,4,2. 
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Table 4.4.1 

RESULTS OF MTCHED-FILTERING 


Prototype : 

C <2> 

Threshold 

Equivalence class obtained 

ICrA 

2, 45, 34, 28 

2(yA 

2, 45, 34, 28, 10, 17, 0, 


20, 41, 5. 

Table 
RESULTS OF 

4.4.'.2 

MATG iED-FILTERING 

Prototype : 

C <3> 

Threshold 

Equivalence class obtained 

la/o 

3, 19, 18, 39, 4^, 44, 45, 


29, 24, 49, 9, 50, 11, 43, 


40, 13.. 



CHAPTER 5 


CONCLUSIONS AInID FURTHER WORK 

Two rather conflicting requirements that are demanded 

of any scheme for representing shape are stability and 

sensitivity, A representation is stable if insignificant, 

similar 

noise-like variations between two/shapes do not cause the 
corresponding representations for them to diverge widely. 
Indeed, stability is essential for the representation to be 
used as a basis on which to classify shapes, in as much as 
it is a measure of correspondence in the representation 
between two similar shapes. On the other hand, the repre- 
sentation should be capable of expressing subtle but 
significant or meaningful differences encountered between 
otherwise similar shapes. In other words, the ease with 
which one can make classifications within a class of 
(grossly) similar shapes is a measure of the ability of the 
representation to capture really fine features of the shape, 
We examined these issues in the context of the Fourier 
representation for shape. We have foun^ that such fine 
distinctions within a class is, indeed, possible by using 
the intuitively satisfying concept of tending-energy, which 
is, related to the Fourier coefficients in a simple way. 



Another useful idea which we explored in this regard is 
that of considering a shape as a signal and employing 
matched-filtering to classify shapes, now considered as 
signals. 

In the context of matched filtering, an idea that 
can be investigated and can prove to be more informative is 
that of the so-called ambiguity function used in radars, 

A plot of the ambiguity function - the ambiguity diagram — 
represents the response of the matched filter to the signal 
for vtiich it is matched as well as to the mismatched signals. 
The ambiguity function is a function of two variables: 
time-delay and frequency. Thus, it is expected that the 
ambiguity function will prove to be more do descriptive 
than the approach we have followed. 

Another question that crops up in judging a .represen- 
tation is: what class of objects are representable? Indeed 
any physically realizable shape is capable of being described 
by a Fourier series representation. One of the objections 
raised against Fourier descriptors is that for shapes with 
’corners’ the Fourier series Converges rather slowly and 
hence the number of significant harmonics one should consider 
to get a specified accuracy of approximation will become 
large. The problem becomes more a ecu te v\hen one considers a 
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class of mixed shapes, some of which are smooth (and hence 
require a lesser number of harmonics) and some of which are 
not so smooth. In such a case, choosing a larger number of 
harmonics might lead to an inefficient solution, from the 
point of view of the time required to calculate the additional 
coefficients, especially if the shapes which have corners 
are few. However, there is a way out; the rate at which 
the coefficients of the Fourier series expansion converge 
to zero for large n, contains information as to the degree 
of smoothness of the curve. 

If a function f(t) and its first k-1 derivatives 

satisfy the Dirichlet conditions and are everywhere continuous 

then as n becomes infinite the coefficients a , and b in the 

n n 

Fourier series of f(t) tend to zero at least as rapidly as 
k+1 

c/n where c is a constant independent of n. If in addition 
the kth derivative of f(t) is not everywhere continuous, then 
either a^ or b^, and in general both, can tend to zero no 
faster than c/n^'*’^. 

Conversely, by observing the rate at which the term in 
the Fourier series of an otherwise unknown function approach 
zero, we can obtain useful information about the degree of 
smoothness of the function, and hence decide whether it would 
be worth our while to consider more harmonics. It is 
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